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We present absorption line analysis of the outflow in the quasar IRAS F04250- 
5718. Far-ultraviolet data from the Cosmic Origins Spectrograph onboard the 
Hubble Space Telescope reveal intrinsic narrow absorption lines from high ion- 
ization ions (e.g., Civ, Nv, and Ovi) as well as low ionization ions (e.g., Cn 
and Si ill). We identify three kinematic components with central velocities rang- 
ing from ~ -50 to ~ -230 km s" 1 . Velocity dependent, non-black saturation is 
evident from the line profiles of the high ionization ions. From the non-detection 
of absorption from a metastable level of C II, we are able to determine that the 
electron number density in the main component of the outflow is ^ 30 cm -3 . 
Photoionization analysis yields an ionization parameter log Uh ~ —1.6 ± 0.2, 
which accounts for changes in the metallicity of the outflow and the shape of 
the incident spectrum. We also consider solutions with two ionization parame- 
ters. If the ionization structure of the outflow is due to photoionization by the 
active galactic nucleus, we determine that the distance to this component from 
the central source is £ 3 kpc. Due to the large distance determined for the 
main kinematic component, we discuss the possibility that this outflow is part of 
a galactic wind. 

Subject headings: seyferts: individual (IRAS F04250-5718); seyferts: absorption 
lines 



Introduction 



Present in the UV spectra of ~50% of Seyfert I galaxies are strong ab sorption troughs 
dDunn et al.l 120071 ) with full- width at half-maximum ~ 20 to 400 km s _1 ( ICrenshaw et al. 



20031 ) that are blue-shifted with respect to their emission line counterparts indicating ma 
terial that is outflowing from the central source. These absorption spectra show lines of 
Civ AA1548,1551; Nv AA1239,1243; Ovi AA1032,1038; and Hi A1216 with lower ionization 



speci es such as Siiv AA1394,1403; and Mgll AA2796,2804 seen less often ( ICrenshaw et al. 



19991 ). Narrow absorption lines (NALs ), defined as h a ying FWHM < 500 km s -1 , are also 
found in quasar spectra. For example, iMisawa et al.l (120071 ) find intrinsic NALs in spectra 
of ~50% of their sample of optically selected quasars at redshi fts z = 2 — 4. Most stud- 
ies of NAL quasars have f o cused on very luminous o bjects (e.g.. iGanguly et al.lll999l . 12001 



Hamann et al.l 1200 ll . 1201 ll ; IMisawa et al.l 120071 120081 ) . NALs provide important diagnostics 
of outflows in active galactic nuclei (AGNs). For example, w e can learn about chemical 



abundances in outflows (e.g., iHamannl Il997l ; iGabel et al.l l2006t lArav et al.l 120071 ) and pro 



vide estimates of their mass flow rates and kinetic luminosities, which are relevent to AGN 
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feedba ck models fe.g..lde Kool et al.ll2001t iHamann et al.ll2001[ iMoe et al.ll2009l ; iDunn et al. 
2010bl lal: Cw et al.lboilh T 



Most likely, the ionization structu re of the outflowing gas is due to photoionization 



rather than collisional ionization (e.g., IHamann fc Ferlandl Il999l ). Therefore, photoioniza 



tion studies are integral to unde rstanding the outflow phenomenon. Computer codes such 



as Cloudy (jFerland et al.lll998l ) have b een developed to self-consiste ntly solve the ioniza- 



tion and thermal equilibrium equations ( jPsterbrock fc Ferlandl 120061 ). Two parameters of 



particular interest are the total hydrogen column density (Njj) of the absorbing material 
and the ionization parameter (Uh), which is proportional to (R 2 n H ) 1 where R is the dis- 
tance to the absorber from the central source, and n H is the total hydrogen number density. 
Thus, determining Uh and measuring n H yield the distance of the outflow. Previous stud- 
ies have found outflows from AGN on scales of ~ 0.1 pc to tens of kpc from the central 



source (e.g.. 


Arav et al. 


Dunn et al. 


2010b). 



Hamann et al. 2001: Kraemer et al. 2001: Moe et al. 2009 



IRAS F042 50-5718 is the first of s ix objects we observed with the Cosmic Origins Spec- 
trograph (COS. lOsterman et al.ll2010[ ) onboard the Hubble Space Telescope (HST), as part 
of our program aiming at determining the cosmological impact of AGN outflows (PID: 11686, 
PI: Arav). IRAS F04250-5718 was chosen for observation because it shows little or no blend- 
ing of absorption features, is bright enough to acquire high S/N in several orbits with the 
high throughput of COS, and has a redshift (z = 0.104; iThomas et al.lll998l ) that allows us 
to cover Lya and Ly/3 as well as C IV, N v, and O VI doublets with the COS far-ultraviolet 
(FUV) g ratings. Observat ion of Lya and Ly/3 is necessary for absolute abundance determi- 
nations ( lArav et al .1120071 ). which will be presented in a future paper. Connected with our 
program, simultaneous XMM-Newton data were acquired (Costantini et al., in preparation). 

IRAS F04250-571 8 is a radio-quiet, low lumi nosity quasar (with absolute blue mag- 



nitude Mb ~ —24.7; IVeron-Cetty &: Veronl 120061 ). Due to the low luminosity and red- 
shift (z = 0.104), IRAS F04250 -5718 is on the quasar /Seyf ert borderline and is spectrally 
classified as a type 1.5 Seyfert (jVeron-Cetty fc Veronl 120061 ) . After some confusion due to 



erroneous positions recorded in some finding charts it was confirmed that the position of 



this galaxy is cons istent with that of the faint blue star LB 1727 (see lHalpern et al.lll998 



Turner et al.lll999l ). Observational history: IRAS F04250-5718 was discovered as an X-ray 
source (1H 0419-57 7) in the First HE AP survey ( [Wood et al.lll984T) and, later, in the Ein- 
stein Slew Survey ( jElvis et al.l Il992l ) and the Swift BAT Survey ( jTueller et al.l 120081 ) . In 
a study identifying optical count erparts to HEAP X-ra y sources, IRAS F04250-5718 was 
identified with a Seyfert galaxy (IBrissenden et al.lll987l ). Pbservations using the Extreme 
Ultraviolet Explorer (EUVE) revealed that this object is extremely bright in the 50-180 A 
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range ([Marshall et al.lll995l h Most of the work on IRAS F04250-5718 to date has focused on 



prop e rties of the X-ray spectral e nergy distribution fe.g.. iGuainazzi et al 



19991 . l2009i ; iPounds et al.ll2004aU bl). Using Suzaku observations, iTurner et al 



1998 



Turner et al. 



(120091 ) found 



evidence of a Compton-thick partial-covering absorber with inferred distances closer to the 
source than the Broad Line Region (BLR). Da ta in the far UV were obtained with th e 
Far Ultraviolet Spectroscopic Explorer (FUSE) ( jDunn et al.ll2007t IWakker fc Savage! 120091 ). 
Optical spectra were obtained using the Australian National University (ANU) telescope 
and the Anglo- Australian Telesco pe (AAT) as well as from the European Southern Obser- 
vatory (ESO) (ITurner et al.lll999l and references therein). To date, no detailed absorption 
line studies in the UV have been reported. 

The plan of the paper is as follows: The observations are discussed in section 2. Column 
density determinations are presented in section 3. In section 4, we discuss photoionization 
analysis of the object, and in section 5, we derive the distance, mass, mass flow rate, and 
kinetic luminosity implied by our analysis. A discussion of our results is given in section 6 
followed by a summary in section 7. The online version of this paper contains a complete 
spectrum of IRAS F04250-5718 with absorption line identifications. 
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2. Observations 

IRAS F04250-5718 was observed for 8 orbits through the COS FUV channel using the 
Primary Science Aperature (PSA) and the medium resolution (R ~ 18000) G130M and 
G160M gratings leading to a total integration time of ~ 1900 s and ~ 15000 s, resp ectively. 
Details about the COS on-orbit performances can be found in lOsterman et al. (bold ). Differ- 



ent central wavelength settings (3 for G130M and 4 for G160M) were used for each exposure 
in order to fill the inter-segment detector gap and minimize fixed-pattern noise, providing a 
continuous spectral coverage from roughly 1130 to 1795 A. 

The reduced data, processed by the COS calibration pipeline CALCOS^, were retrieved 
from the HST archive and combined together using IDL routines^ developed by the COS 
GTO team specifically for COS FUV data. These routines es sentially perform flat- fielding, 



alignment, and co-addition of the processed exposures (see iDanforth et al.l |2010| for de- 



tails). The wavelength calibration of each combined segment is accurate to < 15 km s 1 



( lOsterman et al.ll2010l ). The final combined FUV spectrum of IRAS F04250-5718 rebinned 
to a 20 km s _1 resolution grid has an average signal to noise ratio of ^ 50. Figure [1] shows 
the majority of the COS spectrum (with the original ~ 2 km s~ x oversampling) along with 
identification of intrinsic absorption features. The online version of this paper contains a 
complete spectrum showing details of the line profiles along with identification of Galactic 
ISM absorption features. 

The FUV spectrum of IRAS F04250-5718 displays multiple absorption features from low 
ionization (e.g., Cn, Sim) to high ionization (e.g., Civ, O Vl) ions. Due to the shallowness 
of the Cn trough (residual intensity £ 0.8), we inspected each data set individually and 
found that the absorption feature appears in each subexposure in regions seemingly free 
of flat-fielding issues. We identify three absorption components associated with the UV 
outflow in the un-rebinned COS spectrum of IRAS F04250-5718. These three kinematic 
components, identified using the strong, unblended, Civ and Nv doublet lines, have their 
centroids located at velocities V\ ~ —38 km s _1 , t>2 ~ —156 km s _1 , and v$ ~ —220 km s^ 1 
in the rest frame of the AGN [z = 0.104). 



1 See details in HST Data Handbook for COS at http://www.stsci.edu/hst/cos/documents/handbooks 
/datahandbook/COSJongdhbTOC.html 

2 The routines can be found at http: / / casa.colorado.edu/ ^danforth/costools.html along with an extensive 
description. 
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Fig. 1. — Spectrum of IRAS04250-5718 obtained by COS in June 2010. Absorption troughs 
related to the intrinsic absorber are labeled. Here, ~80% of the spectrum is shown: A full 
identification of all the absorption features is presented in the online version. 
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3. Data Analysis 



3.1. The unabsorbed emission model 



Determining the column densities associated with a given absorption trough requires 
knowledge of the u nabsorbed AGN em ission as well as which emission sources are covered 
by the outflow (e.g., iGabel et al.ll2005l ). The UV emission in AGNs mainly arises from three 



sources; the continuum, the broad emission line region (BELR), and the narrow emission line 
region (NELR) . Afte r correcting the spectrum for Galactic extinction using E(B-V)=0.0134 
( jSchlegel et al.l Il998l ) and Ry = 3.1 and shifting it to the rest frame, the continuum of 
IRAS F04250-5718 is fitted with a single power law of the type F x = Fn 00 (A/1100)^ by 
performing a x 2 minimization over the regions free from known emission or absorption lines. 
A good fit of the continuum of IRAS F04250-5718 was obtained using F 1100 = (3.613 ± 
0.004) x 10~ 14 and = -1.603 ± 0.091. 

In developing an emission model for the BELR and the NELR, we apply the restrictive 
approach outlined in lArav et al.l ( 120021 ) in which the emission lines are modeled using a 
minimum amount of free parameters. We are able to fit the main emission features (Civ, 
O VI, and Lya) using two Gaussian components for the BELs characterized by a full- width at 
half- maximum (FWHM) of ~ 3000 and ~ 9000 km s _1 , though a third, broader component 
(FWHM ~ 14000 km s _1 ) was required to fit the blue wing of the Lya emission. A spline 
fit was used for the Nv line given the difficulty in separating its shallow emission from the 
strong red wing of the Lya BEL. For O VI and C IV, we modeled the NEL of each line of the 
doublet by an additional Gaussian of FWHM ~ 800 km s _1 (fixed by the width of the O VI 
NEL), with the separation of the Gaussians fixed by the difference in wavelengths between 
the doublet lines. A single Gaussian of similar FWHM was used to model the Lya NEL. 
The remaining weaker emission features in the spectrum were modeled by a spline fit. 

The results of the emission modeling process is presented in Figure [2] for C IV and online 
Figures [2jl - [213 for the other strong emission features observed in the IRAS F04250-5718 
spectrum. Individual contributions of each emission source (continuum, BEL, and NEL) as 
well as the total model is plotted as a function of wavelength. In a case where the absorber 
covers only the continuum s ource and the BE LR, we expect to find the absorption troughs 
above the NEL model (see lArav et al.l 120021 ). However, we observe that the flux of the 



deepest Civ absorption trough is lower than each emission component indicating that the 
absorption must cover all three emission sources. 
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Fig. 2. — Details of the unabsorbed emission model for the Civ line in IRAS F04250- 
5718. The total emission model is plotted as a solid line on top of the data. This model 
contains a contribution from three different sources; a continuum (dotted dashed line), a 
broad emission line (dashed line) and two narrow emission lines (dotted lines). The flux of 
the deepest absorption trough is lower than each emission component indicating that the 
absorption must cover all emission sources. 
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3.2. Column Density Determinations 



A first step in characterizing the physical properties of the absorber is the deter- 
minination of ionic column densities associated with the absorption troughs seen in the 
spectrum. Early studies of intrinsic AGN absorption troughs using UV doublets such as 
C IV AA1549,1551 A revealed that they rarely follow the expected 1:2 optical depth ratio be - 
tween the weakest and strongest components of the doublet (e.g. iBarlow fc Sargentlll997l ). 
One interpre tation of this observa t ion is that the absorber does not totally cover the emission 



source (e.g., 



Hamann et al. 



2003: Arav et al. 



2005 



19971 : iGangulv et al.lll99a lArav et aJj [l999| , 120021 : iGabel et al. 



20081 ) . Assuming a partial covering model of a single homogeneous 
emission source, the residual intensity observed for line j of a given ion as a function of the 
radial velocity v can be expressed as: 



/» = /o(«)((l - C ion (v)) + C ion (v)e- T ^) 



X 



where Ci on is the fraction of the emission source covered by an absorber having an optical 
depth Tj, and Iq is the unabsorbed emission. While Equation ([I]) is underconstrained for sin- 
glets (unless assumptions are made about C ion ), multiple lines of a given ion sharing the same 
lower energy level allow us to derive a velocity dependent solution for both C ion and Tj by 
forcing their optical depth ratio R = Ti/Tj to be identical to the one observed in the labora- 
tory R = Xifi/Xjfj. Variations of this technique, including the allowance of di fferent covering 



fractions for each emiss ion source, have been investigated in the literature ([Ganguly et al. 



19991 : IGabel et al.ll2005l ). However, the methods described in those papers cannot be applied 
to doublet lines unless one makes assumptions about the emission source/ absorber model 
(generally assuming the same individual emission source covering for all ions in order to get 
enough constraints on the model). Note that not properly taking into account this partial 



cover ing usually results in an underestimation of the actual column densities (e.g. I Arav et al. 



19991 ). thus hampering the interpretation of the physical properties of the intrinsic absorber. 



While the partial covering model easily accounts for the departure from an optical depth 
ratio of 1:2 between the lines of typical doublets, the physical validity of an absorber model 
displaying a step-function distribution of material accross the emission source is question- 
able. The fact that we typically find different covering fractions for different ions implies 
that the actual distribution of absorber material deviates from this simple absorber model. 
Alternative models in which the emission source is fully covered by a smooth distr i butio n 
of absorbing mat erial accross its spatial extension were studied in Ide Kool et al.l ( 120021 ). 
Arav et al.l ( 120051 ). and lArav et al.l ( 120081 ). Here we use the power law absorber model in 



which the distribution of optical depth accross the source is characterized by t(x) = r max x a . 
In this case, the normalized residual intensity observed for line j of a given ion as function 
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of the radial velocity v can be expressed as 

Ij(v) = [ e~ r — ^ xaimW dx (2) 



Like Equation (JT]), this equation is underconstrained for singlets (unless assumptions are 
made about the shape of the absorbing material distribution parametrized by cij on ) while a 
velocity dependent solution can be derived for both cw and T max ,j for multiple lines of a 
given ion sharing the same lower energy level by forcing their optical depth ratio R = T,i/Tj 
to be identical to the one observed in the laboratory. 

Once an optical depth solution is obtained from either model, the velocity depen dent 



column density is computed using the relation (adapted from lSavage fc Sembachlll99ll ): 



N ion (v) = 3-8 X l ° { Tj (y)) (cm" 2 km" 1 s) (3) 

where Xj is the wavelength of the line (in A), fj is its oscillator strength, and (Tj(v)) is the 
spatially averaged optical depth across the emission source: For the partial covering model, 
( T j( v )) = Ci on (v)Tj(v), and for the power law distribution, (Tj(v)) = T max j(v )/{a ion {y ) + 1). 



3.2.1. Column densities of the metals 

In Figure [3J we show the normalized absorption troughs detected in the IRAS F04250- 
5718 spectrum. The line profiles are presented using the COS sampling of ~ 2 km s _1 . 
In order to compute column densities across the trough, the line profiles are rebinned to a 
common velocity grid with resolution ~ 20 km s _1 , slightly larger than the nominal resolution 
for COS. In determining integrated column densities, we combine kinematic components 2 
and 3 since these components are seen only in Civ and Nv, but not in Lya, Ly/3, or O VI. 
Integrated column densities using the apparent optical depth (AOD) method, the partial 
covering (PC) method, and the Power law model (PL) for components 1 and 2+3 are given 
in Tables [1] and [2J respectively. For doublets, the residual intensity equation (Equation ([I]) 
for partial covering and Equation (j2J) for the power law model) is solved for both lines 
simultaneously. For comparison, we derive column densities using the apparent optical depth 
method (AOD; Ci on = 1) on the weaker line of the doublet. The velocity ranges used for 
integration of column densities are +30 to -90 km s -1 for component 1 and -90 to -290 km s -1 
for component 2+3. The high velocity limit (-290 km s -1 ) is set by the partial observation 
of the Ly/3 line profile in the COS spectrum due to its location at the edge of the detector. 
In §3.2.21 we discuss our treatment of this line. For ions that only appear in component 3, 
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we reduce the integration range to -210 to -290 km s 1 to avoid adding additional noise from 
regions where no trough is detected. 

Absorption lines from Civ AA1548,1551; Nv AA1239,1243; and Ovi AA1032,1038 are 
present in all three kinematic components. Since C IV and N V do not show signs of strong 
saturation, both partial covering and power law methods yield reliable column density mea- 
surements for these ions for each respective absorber model. O VI, on the other hand, shows 
signs of strong saturation not only in the core of each trough, but also in the wings, where 
a 1:1 ratio is observed between the residual intensities of each line of the doublet. Strong 
saturation implies a large deviation of the actual column density in the absorber from AOD 
measurements. The partial covering method also yields only a lower limit to the actual 
column density of O VI due to the upper limit we impose on the optical depth of the weaker 
line: We limit the optical depth of the weaker line to r < 4 because the difference between 
residual intensities of the lines is well within the noise whenever the optical depth of the 
weaker line of the doublet is ^ 4. Similar constraints are imposed on the power law method 
((t) < 4 for the weaker line of the doublet) resulting in a column density that is only ~20% 
larger than that estimated by the partial covering method. 

Absorption lines from Cll A1335; Sim A1207; Siiv AA1403,1394; and S IV A1063 are 
present only in kinematic component 3. For Si IV, we find that the residual intensity of the 
blue line of the doublet closely matches that of the red line when properly scaled by the ratio 
of line strengths (Is = Ir), indicating that the AOD method yields reliable optical depths 
for this shallow trough. We then assume that this is also the case for the even shallower 
troughs from C II, Si ill, and S IV. 
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Fig. 3. — Normalized absorption troughs detected in the IRAS F04250-5718 spectrum. The 
strongest line of a doublet/multiplet is plotted in black while the weakest is plotted in gray 
(red in the online version). For each doublet line as well as the Hi lines, we indicate the 
range of physically allowed residual intensity for the strongest line by plotting a dotted line 
showing the expected residual intensity for the strongest component if the emission source 
was fully covered by the absorber (optical depth ratio of R = A»/»/ Xjfj)- The Ly/3 is only 
partially detected in our COS observations due to its location at the edge of the detector. 
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3.2.2. Column density of Hi 

The Ly/3 trough is located near the edge of the blue detector (~1131 A) and only 
partially observed. Given the poor wavelength solution at these low wavelengths, we shifted 
the entire Ly/3 trough by 39 km s _1 to match the kinematic structure of Lya. The poor 
wavelength solution may affect the shape of the line profile as well. In order to assess 
whether the Ly/3 trough can be used in our analysis, we compared the COS data with 
2006 FUSE observations of IRAS F04250-5718. The FUSE spectra were downloaded from 
the Multimission A rchive at Space Telescope (MAST) and processed with CalFUSE v3.2.3 



( iDixon et al.l 120071 ) . Due to low photon counts, we opted to use bot h the night and day 



exposures in the coadded spectrum since the dayglow emission features (IFeldman et al.ll200ll ) 
did not affect our scientific goals. We flux calibrated all eight segments to match the flux 
level in the LiF la segment. Comparison of the overlapping regions of the 2010 June COS 
and 2006 May FUSE spectra shows no variations in the continuum level or line profiles at 
the S/N of the FUSE data. Similarly, the shapes of the Ly/3 troughs in the two instruments 
are identical within the noise. We conclude that a simple wavelength shift of the Ly/3 trough 
in the COS data suffices for our analysis. 

Hi column densities are derived by treating the Lya and Ly/3 lines as a doublet, and 
comparisons are made to AOD measurements using Ly/3. The lines are clearly separated 
in component 1 allowing for reliable estimates of the H I column density in this component 
for respective absorber models. In component 2+3, the deepest parts of the Lyct and Ly/3 
troughs look close to saturation, having a nearly 1:1 residual intensity ratio. Inspection of 
the Lyman series lines in the FUSE spectrum indicates a possible saturation of the Ly/3 
and Lyy lines. Further, it is difficult to estimate the unabsorbed emission in the region of 
the Ly/3 absorption since only part of the emission is observed and the COS S/N is severely 
degraded at the edge of the detector. We therefore consider the column density of H I derived 
from the partial covering method a lower limit for kinematic component 2+3. The column 
density derived from the power law method is ~ 50% higher. The lower S/N of the FUSE 
data does not allow us to get a better estimate of the Hi column density or determine a 
better unabsorbed emission model. However, the absence of a bound-free edge for H I in the 
FUSE data allows us to place an upper limit on the Hi column density of 10 15 ' 8 cm -2 . 
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Table 1. Column densities for component 1 of IRAS F04250-5718 



Ion N ion (AOD) N ion (PC) N ion (PL) 



Hi U2±H 532±|g 1340±f° 

Civ 79.5±i'| 113±^ 3041?, 06 

Nv 59.7t\i 73.2ti° 7 4 135lg 

Ovi 63411? 1220±ao 2660± 090 



Notes. Column densities are given as linear values in units of 10 cm . Asym- 
metrical er rors are from phot on statistics only and are derived using the technique 
outlined in iGabel et al. ( 2005 ). 



Table 2. Column densities for component 2+3 of IRAS F04250-5718 



Ion N lon (AOD) N lon (PC) N lon (PL) 



Hi 


17301™ 


21001^° 


3 070 +990 


Cii 


n-i±i:i 






Civ 


472 ±\ 


610lf 


1160lf° 


Nv 


637±1 


834+^ 


1360lf° 


Ovi 


27OOI20 


3390l 2 ° 40 


417Q +1330 


Si in 


9 44+ - 06 






Si iv 


11 q +0 - 7 


J -*- L -0.6 


i0 - l -o.6 


Siv 


101II 







Notes. Column densities are given as linear values in units of 10 12 cm~ 2 . Asym- 
metrical er rors are from phot on statistics only and are derived using the technique 
outlined in IGabel et al.l ( 2005 ). 



- 15 - 



Fig. 4.— Comparison between the FUSE and HST/COS spectrum of IRAS-F04250-5718 
taken four years apart. No noticeable changes in the continuum level or absorption line 
profiles is observed within the FUSE S/N between these two epochs though some variations 
may have occured between observations. The Ly/3 trough observed in the FUSE data has 
a line profile that is identical within the noise to the one observed in the COS observation, 
suggesting that the COS wavelength solution in this region can be, as a first approximation, 
corrected by a simple shift of the line profile. 
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3.2.3. Effects of Deconvolution on Column Density Determinations 

Early on-orbit analysis of the COS Line Spread Function (LSF) revealed a significant 
departure from the Gaussian profile observed during ground testing, showing a broadened 
central core associated with strong broad wings scattering a significant part of the flux 
(e.g. Kriss 2011 COS Instrument Science Report 2011-01(vl))§. The major effect of such a 
degraded LSF is to scatter photons from the continuum region into the absorption troughs, 
decreasing their apparent depth and hence affecting the determined column density derived 
for each ion, especially for narrower absorption troughs with Av ~50 km s _1 (Kriss et al. 
2011, accepted for publication in A&A). 

Given the high signal to noise of our spectrum, we attempted to correct for the degraded 
LSF by deconvolving the spectrum using the Lucy- Richardson algorithm implemented in the 
STSDAS routine lucy in the analysis . restore package. Given the wavelength dependance 
of the LSF, the spectrum was deconvolved in 50 A intervals around the observed absorption 
troughs using the appropriate LSF for each region (see Kriss et al. 2011 for details). Column 
densities are then computed for the deconvolved absorption troughs after rebinning the data 
to the identical 20 km s _1 velocity grid. For the partial covering method, the typical increase 
in the computed column density for the broadest troughs (e.g. component 2 of Civ) is of 
the order of several percent, but can reach up to ~30% for the narrower troughs such as 
Si IV. While these results show a general increase in the measured column densities on the 
deconvolved line profile, the significant increase of the associated noise in the continuum 
implied by the deconvolution process and the inability to accurately characterize that noise 
increase lead us to consider the column densities derived on the non-deconvolved spectrum 
in the remaining analysis while bearing in mind that the column densities associated with 
narrow troughs can be underestimated by up to ~ 30% for the partial covering solution due 
to the degraded LSF. 



3 found at http:/ /www.stsci.edu/hst/cos/documents/isrs/ISR2011_01(vl).pdf 
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4. Photoionization Analysis 

4.1. Spectral Energy Distribution 

The photoionization and thermal structure of an outflow depends on the spectral energy 
distribution (SED) incident upon the outflow. A st andard SEP (hereaf t er MF 87) used in 



photoionization modeling of AGN was developed in iMathews fc Ferlandl (119871 ). This SED 



contains a prominent "big blue bump," which came under scrutiny in later works. For 



example, lLaor et all (119971 ) dispute the existence of an X-ray-to-ultraviolet (XUV) bump in 
AGN based on ROSAT data. H owever, X-ray data of IRAS F04250-5718 seem to require an 
XUV bump jTurner et al1ll999h . Since this object's SED is unobservable in the UV ionizing 



region, we analyze models with different XUV bumps. Studies of the SED in the X-ray 
regime have shown extreme variability (e.g. Guainazzi et al. 1998). Pounds et al. (2004a) 
showed that the variability could be explained by absorption of a relatively constant SED by 
cold dense matter. In order to ameliorate this issue, we use the simultaneous observations 
from HST/COS and XMM-Newton to constrain the SED of IRAS F04250-5718: HST/COS 
data constrain the SED in the frequency range 2 x 10 15 Hz £ v £ 3 x 10 15 Hz, while 
data from XMM-Newton (Costantini et al. 2011, in preparation) constrain the high energy 
portion (y ^ 10 17 Hz) of the SED incident upon the UV absorber. While the FUSE data 
were taken four years earlier, they agree with the COS data in the overlap region. The 
co-added COS and FUSE data can be fit with a single power law (see section 3.1) extending 
the far UV constraints to ~1 Ry. 

We consider two families of SEDs in the analysis of IRAS F04250-5718 for a total of 
three SEDs. Figure [5] shows the three SEDs along with the MF87 SED for comparison. For 
two of our SEDs, we use piecewise power law models. The blue line in Figure [5] (hereafter 
BPL1) is the softest SED we consider physically reasonable given the constraints imposed 
by the data. It has a spectral break at ~1 Ry, very near the edge of detection in the far UV. 
The red line in Figure [5] (hereafter BPL2) represents the broken power law model that best 
fits the data while containing an blue-bump that is no more prominent than that in MF87. 

For our final SED (hereafter AGN-UVbump), we generate the sum of two continua using 
the agn command in Cloudy. This command produces the continuum given by 

f v = az/^exp f ~ V ) + bu a % (4) 



kT, 



BBS 



where v is the frequency, k is Boltzmann's constant, and h is Planck's constant. The spectral 
index a uv gives the low-energy slope of the optical/UV component, a x gives the X-ray 
slope, and Tbbb is the temperature of the "big blue bump." The coefficient mutliplying 
the optical/UV power law component provides an exponential cutoff for energies ^ 0.01 Ry, 
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and the coefficient multiplying the X-ray power law component is adjusted to produce the 
optical to X-ray slope (a ox ) provided by the user. The X-ray component is set to zero for 
energies < 0.1 Ry and falls off as v~ 2 for energies £ 100 keV. Generated by the Cloudy 
command agn 2.8e6 -1.30 -0.80 -0.63 (where the numbers are Tr RR , a nx , a ,,,, ., and a x , 
respectively), the SED AGN-UVbump is similar to the one used by lArav et all ( 120071 ) for 
Markarian 279 with the main differences being the frequency at which the spectral break of 
the big blue bump occurs a nd the hard X-ray excess seen in the XMM-Newton data (also 
observed in Suzaku data by lTurner et al.l 120091 ) . 



Table |3] shows the bolometric luminosity (L^i) and the number of hydrogen ionizing 
photons emitted by the source per second (Qh) for each SED. 



Table 3. SED Emission Properties 



SED L bol (10 45 erg/s) Q H (10 55 s" 1 ) 



MF87 6.4 5.2 

AGN-UVbump 8.8 5.0 

BPL1 6.6 3.6 

BPL2 7.3 5.0 
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Fig. 5. — SEDs used in the analysis of IRAS F04250-5718. Ionization potentials of Hi and 
O VI are marked with dotted black lines to mark the UV ionizing region of the SEDs. The 
MF87 SED drops to a minimum of about 10 43 5 erg s" 1 at E m 365 eV. 
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4.2. Photoionization Modeling 

The total hydrogen column density Nh and ionization parameter 

47rit z ri H c 

for the absorber are determined by self-consistent photoionization modeling (where c is the 
speed of light, Q H is the rate of hydrogen ionizing photons emitted by the AGN, R is the 
distance to the absorber from the central source, and nn is the total hydrogen number den- 
sit y). For this, w e use v ersion c08.00 of the spectral synthesis code Cloudy, last described 



by iFerland et al.l (119981 ) . In these models, we assume solar abundances as given in Cloudy 
and a plane-parallel geometry for a gas of constant hydrogen number density. In §4.2.1[ we 
discuss the effects of supersolar metallicities on our results. In determining the ionization 
parameter that leads to our distance determination, we assume that the high ionization fea- 
tures (Civ, Nv, Ovi) come from the same gas as the lower ionization species (Cll, Si ill, 
Si IV, S iv). This is supported by the kinematic correspondence for all of the ions in detected 
in the outflow (see Figure [3]). 

To determine Nh and Uh, we g enerate grids of m odels varying these parameters in 0.1 



dex steps (similar to the approach of lArav et al.ll200ll ). For a given SED, ~4500 models are 



computed to investigate a parameter space with 15 < log Nh < 24.5 and -5 < log Uh < 
2. At each point in a grid, we tabulate the predicted column densities of all relavent ions 
and then compare the measured column densities to those predicted by the models. Figures 
[6] and [7] show results of this analysis for kinematic components 1 and 2+3 of the outflow, 
respectively, using the SED AGN-UVbump developed in §4.11 In these figures, lines where 
a linear interpolation of the tabulated Ni on matches the measured Ni on are plotted in the 
N h -Uh plane. 

In spectra with absorption features from three or more lines from t he same energy level 



of the same ion, it is possible to test the different absorber models (jde Kool et al. 



2002 



Arav et al.l 120051 ). However, in IRAS F04250-5718, we have no observational constraints 
on the spatial distribution of absorbing material across the emission source. Therefore, we 
use the average value of N ion determined by partial covering and power law methods for 
Civ, Nv, ,Ovi, Hi, and SilV. For the remaining ions (Cll, Si ill, and Siv), we have ionic 
column densities determined by the AOD method only. The fact that all three measurement 
methods give similar results for the column density of Si IV indicates the AOD method yields 
a reliable estimate of the actual column density for this ion. We assume this is the case for 
the other ions with shallower troughs since the line profiles are similar. 



The first kinematic component has only four constraints since none of the low-ionization 
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species are detected there. An approximate solution is determined by visual inspection of 
the graph and plotted as a solid square in Figure [6] at log Uh ~ -1-4 and log N H ~ 18.9. 
The error bars span the values of Njj and [/# that yield N ion to within a factor of ~ 2. The 
parameters are correlated such that increasing Uh requires an increase in Nh resulting in a 
line of solutions that runs roughly parallel to the curve plotted for Nv from the center of 
the solid square. Table H] details the results of models using different SEDs for component 1 
of the outflow. 

Figure [7] shows the results for kinematic component 2+3 using the SED AGN-UVbump. 
Four different solutions are given: The solution marked with a plus is approximately equidis- 
tant from the C II, C IV, and O VI curves and encompasses all lines to within a factor of ~ 2.5. 
The crossing point of Si ill and Si IV is marked with a square, and the crossing point of C II 
and C IV is marked with a diamond. Finally, our best by-eye fit to the data is marked with 
a triangle at log Uh ~ —1.8 and log Nh ~ 19.6. The "equidistant fit" and the "best by-eye 
fit" both provide reasonably good fits to the data. We prefer the best by-eye fit since it 
obeys the firm lower limit determined by iV(0 vi) while fitting all of the ions except C II to 
within a factor of ~ 2.2 as well as improving the fits to Nv and Siiv, while giving results 
for Civ, Si hi, Siv, and Hi that are similar to those given for the equidistant fit. The 
disadvantage of the best by-eye fit is that iV(C n) is underpredicted by a factor of ~ 5 while 
for the equidistant fit, iV(Cll) is underpredicted by a factor of ~ 2.2. Table |5] details the 
results of models using different SEDs for component 2+3 of the outflow. 

Solutions for all three SEDs are compared in Figure [8J A dashed box encloses the 
solution space that gives reasonable fits to the data while allowing for changes in SED shape 
and chemical abundances. Solutions determined by the crossing of C II and C IV in the iV#- 
Uh plane are omitted because they severely underpredict O VI, for which we have a firm lower 
limit. Less severe underpredictions may be alleviated by supersolar metallicities, which we 
investigate in the next section. 
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Fig. 6. — Component 1: Plots of constant N ion in the Nh-Uh plane. The plotted lines 
trace models whose predicted N ion matches the observed values. The solid black square 
marks the model that best fits all of the given constraints. The horizontal dashed line at log 
Nh = 24 is the boundary above which an absorber becomes optically thick due to electron 
scattering. The slanted dashed line is the approximate location of the hydrogen ionization 
front. The error bars mark solutions that are within a factor of ~ 2 for all ions. Nh and Uh 
are correlated such that increasing Uh requires an increase in Nh with the line of solutions 
running roughly parallel to the Nv line from the center of the filled square. 
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Fig. 7. — Component 2+3: Plots of constant N ion in the Nh-Uh plane. The plotted lines 
trace models whose predicted N ion matches the observed values. Symbols mark the models 
that fit certain constraints. The plus is equidistant from C II, C IV, and O VI. The square 
and diamond mark the crossing points of silicon ions and carbon ions, respectively. The best 
by-eye fit is marked with a triangle. The lower limit on H I is due to the fact that Ly/3 may 
be saturated, and the upper limit is determined from non-detection of a Lyman break (see 
§3.2.2p . The shaded band therefore represents values of Hi consistent with the data. 
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Table 4. Photoionization models for component 1 (-90 km s 1 < v 

< 30 km s" 1 ) 



Ion log(V lo „) (cm- 2 ) lo §(fe) 

observed t AGN-UVbump BPL1 BPL2 



log U H ■■■ -1-40 -1.23 -1.21 

log N H ■■■ 18.85 18.90 18.87 



H I 


14.97j££ 


+0.02 


-0.21 


-0.19 


C IV 


14.32t°i? 


-0.21 


-0.21 


-0.21 


N V 


i4.02±g:|| 


+0.18 


+0.20 


+0.20 


VI 


15.29_ ; 20 


-0.20 


-0.17 


-0.19 



' Average of column densities derived from the partial covering and power law meth- 
ods. Uncertainties reflect the range of column densities with the higher given by the 
power law method and the lower given by the partial covering method. 
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Table 5. 



Photoionization models for component 2+3 (-290 km s 1 < v 
< -90 km s" 1 ) 



Ion 


log(V io „) * (cm- 2 ) 
observed 


AGN-UVbump BPL1 


BPL2 


log U H 
log N H 




-1.75 -1.60 
19.58 19.59 


-1.57 
19.57 



H I * 


15.41-15.82 


+0.31 


+0.06 


+0.08 


C II * 


lO.UO_ 05 


-0.42 


-0.50 


-0.67 


C IV t 


14.951°;^ 


+0.29 


+0.29 


+0.28 


N V t 


15.041°;°! 


+0.02 


+0.01 


+0.01 


VI t 


1 5 ^8+ - 04 

lO.OO_ 05 


+0.01 


-0.01 


-0.03 


Si III * 


12 39+ 001 


+0.03 


+0.04 


+0.03 


Si IV t 


1Q 1 fi +0.02 


+0.05 


-0.03 


+0.06 


S IV * 


14 00+ 002 

i ^' uu -0 .01 


-0.10 


-0.14 


-0.11 



* Upper limit determined from the non-detection of a Lyman break in FUSE data. 
The range of values for H I is shown as a gray strip in Figure [7] 

* Column densities derived from the AOD method. Uncertainties are due to photon 
statistics only and are thereby underestimated. 

' Average of column densities derived from the partial covering and power law meth- 
ods. Uncertainties reflect the range of column densities with the higher given by the 
power law method and the lower given by the partial covering method. 
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4-2.1. Supersolar Metallicity Models 

The ionization parameters for the outflow were determined under the assumption of solar 
abundances. The grids of models support this assumption for kinematic component 1 of the 
outflow since the metal lines have crossing points very close to the Hi line, but kinematic 
compononent 2+3 shows signs of supersolar metallicities. In particular, the equidistant 
solution and the best by-eye fit solution both overpredict H I: In Figure [7J these solutions 
lay above the line marking the Hi upper limit. In order to investigate the effects of changes 
in metallicity (Z) o n our determination of Nh and Uh, we use the abundances given by 



Ballero et al.l ( 120081 ) for Z/Z Q =4.23 and 6.11. They also list abundances for Z/Z Q =7.22, 
but these do not yield a good solution for any of our SEDs. The results are plotted in 
Figure [8l For each SED, there is an increase in Uh of ~0.1 dex and a decrease in Nh of up 
to ~0.4 dex. We note that the solutions for Z/Zq—A.23 were generally better than those 
for Z/Zq— 6.11, which underpredicted Hi for both of the broken power law SEDs. For the 
AGN-UVbump SED, the solution found using Z/Z & =6.11 is plotted but left out of the box 
of reasonable solutions since the fit was so poor. The C IV was fit better in the supersolar 
models than in the solar ones due to a relative increase in abundance of silicon and nitrogen 
with respect to carbon, but this resulted in an underprediction of C II that is more severe 
than in the solar models. In fact, for the AGN-UVbump SED, using Z/Z Q =4.23, all lines 
are fit to within a factor of ~ 1.1 except for Cn whose column density is underpredicted by 
a factor of ~ 10. We defer the full determination of chemical abundances for the outflow in 
IRAS F04250-5718 to a future paper. 
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Fig. 8. — Photoionization solutions for component 2+3 of the outflow for different SEDs and 
metallicities. The colors match those of the SEDs in Figure [5j Green is the AGN-UVbump 
SED, blue is BPL1, and red is BPL2. Abundances are solar unless noted otherwise: Solar 



abundances are from Cloudy, and supersolar abundances are from iBallero et al.l ( 120081 ). 
Only the best by-eye fits are shown for supersolar metalllicites. The dashed box shows the 
range of plausible solutions. 
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J^.2.2. Multiple Ionization Parameter Models 

For optically thin models, multiple ions from the same chemical species offer (ap- 
proximately) metallicity independent solutions for Nh and Uh- In the COS spectrum of 
IRAS F04250-5718, we have C II, C IV, Si ill, and Si IV. The results obtained using just the 
ratios of N(C il) / N(C iv) and iV(Si ill) / iV(Si iv) are shown in Figure El The diamonds mark 
the results from carbon and the squares mark the results from silicon. While the values of 
Nh agree very well, there is a discrepancy in the determination of Uh of ~ 0.4 dex. Fur- 
thermore, the ratio N(Ctv)/N(Cu) is substantially overpredicted by the models that give 
reasonable fits to the other ions. We therefore investigate the possibility of two ionization 
parameters improving on the photoionization solutions for kinematic component 2+3. 

For all three SEDs, the by-eye best fits (marked with triangles in Figure |8]) underpredict 
the column densities of Cn and Sim. By choosing a second ionization parameter that is 
near the crossing point (log Ujj ~ —3.5) of Cn and Sim in the N h -Uh plane, we improve 
the fits to these two low ionization ions while leaving essentially unchanged the predicted 
column densities of the other ions when values from the two models are summed. However, 
low ionization parameter models predict a substantial amount of Si II, which we do not detect 
in the data. Using an upper limit of log iV(Sill) ~ 11.9 cm -2 , the two ionization parameter 
solutions overpredict iV(Sill) by ~ 0.6 dex. We conclude that the two ionization parameter 
does not significantly improve our results, simply trading the underprediction of N(C il) for 
the overprediction of iV(Sill). 
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5. Absorber Distance and Energetics 

The distance to the absorber can be computed from the ionization parameter (see Equa- 
tion ()5])) if the total hydrogen number density n H is known. To estimate n H in kinematic 
component 3 (the high velocity component) we consider the ratio of populations of the first 
excited metastable and ground levels of C II. The rebinned line profiles of C II A1335 and C II* 
A1336 are shown in Figure!! Integrating the flux in the region where a C II* trough should be 
according to the detected C II trough yields a column density for C II* of 2.81^ 1 ' 7 x 10 12 cm -2 
where the errors are 1-cr photon statistics. Due to non-detection of the excited state of C II, 
we place a conservative upper limit on the column density of Cll* of 6.74 x 10 12 cm -2 by 
assuming the noise could hide up to a 3-cr detection. Comparing collisional excitation models 
with the ratio of the column density measurement of Cll (11.1 ± 1.3 x 10 12 cm -2 ) and the 
upper limit on the column density of Cll* ( £ 6.74 x 10 12 cm -2 ), we place an upper limit on 
the free electron number density n e ^ 30 cm -3 (see Figure [TO]) . Since our photoionization 
modeling implies we are within the He ill region, n e « 1.2n H . In order to place robust lower 
limits on the distance, mass flow rate, and kinetic luminosity of the absorber, we adopt the 
model with the SED BPL2 and a metallicity of Z = 4.23 which yields the maximum Uh 
(10 -14 ) and minimum Njj (10 19 3 ) of all our models. These results imply a distance of 

R = ( -. Q T H T ) ' £ 3 kpc (6) 

for kinematic component 2+3. We note that the data hint at a detection of Cll*. Since the 
Cll* trough that appears in Figure is too shallow to clearly distinguish from the noise, we 
do not claim we have a measurement of the column density of C II*, but state that the lower 
limit on the distance may be an estimation of the actual distance. We cannot determine the 
distance to kinematic component 1, since we have no density diagnostics for this component. 

In order to estimate the mass, mass flow rate, and kinetic luminosity of the outflow, we 



must assume a geometry. For a full discussion of the geometry used here, see lArav et al. 



( 120111 ). Assuming that the outflow has the geometry of a thin, partially-filled shell, the mass 
is given by 

M out = 4nfim p R 2 N H tt > 2 x lO 7 fi o . 5 M , (7) 

where \x = 1.4 is the mean atomic mass per proton, m p is the mass of the proton, and 
Q is the covering fraction of the outflow as seen from the central source. We adopt the 
value Q = 0.5 due to the fact that outflows are seen in about 50% of observed Seyfert 



galaxies (IDunn et al.l 120071 ) and the luminosity of IRAS F04250-5718 is just above that of a 



Se yfert galaxy. A com prehensive discussion of covering fraction for quasar outflows is given 



in 



Dunn et al.l (l2010bl ). In order to reflect the uncertainty in Q, we provide our results in 
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units of = fi/0.5. The average mass flow rate is just the total mass divided by the 
dynamical timescale R/v: 

M out = 4irfim p RN H vn > 1 Q . 5 M^r" 1 . (8) 
This value is similar to the mass accretion rat e we find unde r the assumption that the black 



hole is accreting mass at the Eddington rate (jSalpeterl Il964l ) : 



Lboi 



M acc = ~ 2 Moyr" 1 , (9) 
ecr 

where we have assumed an efficiency e = O.f. Finally, the kinetic luminosity is given by 

E k = ^f- > 1 x 10 40 erg s~K (10) 

In Table El we compare these results with those we have found in high luminosity quasars. 

Given that the bolometric luminosity is Lj, i ~ 9 x 10 45 erg s _1 , this outflow is not 
energetically signific ant in AGN feedback scen arios, which typically require Ek to be a few 



percent of L^oi (e.g., IScannapieco &: Ohll2004f ). The mass flow rate is also rather small for 



the purpose of chemically enriching the environment. Assuming a duty cycle of 10 8 for the 
active phase of this AGN, Hellman and Arav (2011, in preparation) find that a mass flow 
rate of 1 solar mass per year is only enough to enrich a small galactic halo of ~ 1O 1O M to the 
observed Intra-Cluster values, or alternatively 1 Mpc 3 of the Inter-Galactic Medium. The 
outflows observed in the high luminosity objects shown in Table [6] are much more important 
for AGN feedback, both energetically and for chemical enrichment. 



Table 6. Properties of Large Scale Outflows Measured by Our Group to Date 



Object 


log L bol 
(erg s" 1 ) 


log ra e 
(cm~ 3 ) 


v t 
(km s _1 ) 


Av * 
(km s — 


IRAS F04250-5718 


45.9 


< 1.4 


-200 


200 


QSO 1044+3656 


46.8 


3.8 


-3995 


420 


AKARI J1757+5907 


47.6 


< 3.8 


-925 


250 


SDSS J0318-0600 


47.7 


3.3 


-4160 


890 


SDSS J0838+2955 


47.5 


3.8 


5000 


2000 


QSO 2359-1241 


47.7 


4.4 


1385 


130 



R log Nh log Uh log -Efc M ou t References 

(kpc) (cm -2 ) (erg s _1 ) (Mq yr" 1 ) 



> 3.7 > 19.3 < -1.4 > 40.13 > 1.07 1 
1.7 ± 0.4 20.8 ± 0.1 -2.19 ± 0.10 44.81±[5 i? 120 ±25 2 

> 3.7 > 20.8 > -2.2 > 43.30 > 70 3 
5.9 ± 0.4 19.9 ± 0.2 -3.1 ± 0.1 44.551°; .H 60 ± 20 4 



3.3±{ l 20.8 ± 0.3 -2.0 ± 0.2 45.351°; H 3001^0 
3.2±^;f 20.6 ± 0.2 -2.4 ± 0.2 43.36 ± 0.27 90I20 



"iMean velocity of the kinematic components used in the analysis, 
t Full-width at zero-intensity. 

Arav et alj [120111) , (3) lAoki et all d201lh . (4) bunn et"al] i bOlObl '). (5l lMoe et ail j2009h . (6l lArav etaLI 



*(1) This work, (2 
(2008): [Ko"rista et alj fcoosh 
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Fig. 9. — Line profiles for Cll and Cn* rebinned to 20 km/s. We place an upper limit on 
C II* by assuming the noise could hide up to a 3-a detection. The dashed line shows the C II 
trough scaled to the 3-a noise in the region where C II* would be. 
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lo gio( n e) ( cm 3 ) 

Fig. 10. — The ratio of the level population of the first excited state (E = 63 cm *) to 
the level population of the ground state for Cll versus electron number density (n e ) for T 
= 10000 K (dotted line), T = 15000 K (solid line), and T = 30000 K (dashed line). Non- 
detection of Cll* implies a level population ratio of ~ 0.25 ± 0.12. For a conservative upper 
limit, we use the 3-a value of 0.61 yielding n e < 30 cm -3 . 
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Discussion 



6.1. Comparison with Galactic Winds 



While it is clear that outflows at distances of 



< 



10 pc and/or with velocities ^ 500 km s 1 
(jMartinl 120051 ) are AGN driven, outflows at distances £ 1 kpc with lower velocities, such 



as the one we find in IRAS F04250-5718, are difficult to distinguish from galactic winds. 
Galactic winds are ubiquitous in starbur st galaxies: Studies show ^ 75% of ultra-luminous 
infrared galaxies (ULIRGs) have winds ( Veilleux et al. 20051 an d references therein). An op- 
tical study of edge-on Seyfert galaxies by I Colbert et al.l (119961 ) revealed large-scale outflows 
in ~25% of these AGN. However, since evidence of starburst activity is fo und in about half 
of optically selected Seyfert 2 galaxies ( iGonzalez Delgado fc Perezl |2000| ). it is difficult to 
determine whether the galactic-scale outflows in Seyferts are driven by the AGN or star- 
bursts. In a sample of ultra-lu minous infra-red gala xies (ULIRGs) containing both an AGN 
and a high star formation rate, iRupke et al.l (120051 ) found that the velocity distributions in 
outflows observed in Seyfert 2 galaxies differed from those found in starburst galaxies with 
no observable AGN activity, suggesting the AGNs play a role in driving the galactic-scale 
outflows in these objects. 

Much of the work on galactic winds has focused on optical em ission lines produced 



by the warm gas (e.g., IColbert et al. 



1996 



Lehnert fc Heckmanlll996l ). and X-ray emission 



lines produced by the hot gas (e.g., iDahlem et al.lll998l ). Emission line techniques favor 
edge-on galaxies where it is easier to separate the wi nd emission from background emis 



sion. Complementary s tudies of absorption lines (e.g., iHeckman et al.ll2000l ; iMartinl 12005 
Martin &: Bouchel |2009| ) favor face-on galaxies. Absorption line studies of galactic winds 
afford the advantage that the whole range of gas densities in outflows can be pro bed unlike 
emission line studies where the densest gases are favored (IHeckman et al.ll2000l ). Much of 



this wor k has relied on N a l AA5 892,5898 lines and find they do not completely cover the 



source. iMartin fc Bouchel (120091 ) found they could better constrain the kinematics of the 
outflow using Mgll AA2796,2803 since the absorption features are well separated allowing 
for better analysis than the heavily blended Nal doublet. 

T he distance and spe ed we find for kinematic component 2+3 is typical for a galactic 
wind (jVeilleux et al.l 120051 ) . We are unable to determine whether this outflow is AGN or 
starburst driven, but partial-covering of the emission sources implies the outflow is intrinsic 
to the host galaxy so photoionization of the outflow is likely dominated by AGN emission. 
In this case, regardless of outflow's origin, photoionization analysis, along with the upper 
limit on the electron number density provided by a detection of Cll with no detection of 
Cll*, yields a lower limit on distance to the absorber from the central source. 
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7. Summary 

The HST/COS spectrum of the low luminosity quasar IRAS F04250-5718 reveals an 
outflow with both high and low ionization species. The lines are narrow (FWHM £ 200 km 
s _1 ), and similar redshifts for corresponding emission and absorption features implies the 
outlflow is intrinsic to the AGN. Further evidence of the intrinsic nature of this outflow is 
partial covering of the emission source evident in Civ, Nv, Ovi, and Hi. 

Using the measured column density of Cll along with an upper limit on the column 
density of Cll*, we found a conservative upper limit on the electron number density in the 
outflow of ~ 30 cm~ 3 . Photoionization modeling yielded a range in ionization parameter 
log Uh ~ — 1.8 ± 0.2 and total hydrogen column density log N H ~ 19.55 ± 0.25 cm -2 . Using 
model assumptions that minimize the distance to the outflow, we determined a conservative 
lower limit of R ^ 3 kpc. Assuming a thin-shell geometry, lower limits on the mass flow rate 
and kinetic luminosity were found to be ^ 1 M Q yr _1 and ^ 1 x 10 40 erg s _1 , respectively. 

In determining limits on the distance, mass flow rate, and kinetic luminosity of the 
outflow, we investigated the effects on our results of changing the SED incident upon the 
absorber and the metallicity of the gas as well as differences in column density determinations 
when different absorber models (i.e. partial-covering versus power-law models) were assumed. 
The ionization parameters found for the three different SEDs are within ~ 0.2 dex of each 
other, while increasing the metallicity resulted in ionization parameters ~ 0.2 dex higher 
for each SED and lowered Nh by ~ 0.3 dex from the best-fit solar models. Changing the 
SED had little effect on the total hydrogen column density, while differences in the absorber 
model had little effect on ionization parameter and changed N H by < 0.1 dex. 

We acknowledge support from NASA STScI grants GO 11686 and GO 12022 as well as 
NSF grant AST 0837880. We would also like to thank Pat Hall for insightful suggestions 
and discussions. JIGS and CB acknowledge financial support from the Spanish Ministerio 
de Ciencia e Innovacion under project AYA2008-06311-C02-02. 

when changes in SED shape and chemical abundances are considered as well as differ- 
ences in column density determinations when different absorber models (i.e. partial-covering 
versus power-law models) are used 
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